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Abstract —The statistical properties of the multivariate Gamma- 
—Gamma (FF) distribution with arbitrary correlation have remained 
unknown. In this paper, we provide analytical expressions for the 
^ joint probability density function (PDF), cumulative distribution 
^-^function (CDF) and moment generation function of the multivariate 
FF distribution with arbitrary correlation. Furthermore, we present 
novel approximating expressions for the PDF and CDF of the sum 
of FF random variables with arbitrary correlation. Based on this 
statistical analysis, we investigate the performance of radio frequency 
and optical wireless communication systems. It is noteworthy that 
the presented expressions include several previous results in the 
I— iliterature as special cases. 

[- I Index Terms —Diversity receivers, free-space optical communica¬ 
tions, multivariate FF distribution. 
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I. Introduction 

In practical wireless communications, the electromagnetic sig¬ 
nal experiences composite small-scale fading and shadowing 
simultaneously. Over the past decades, several works have fo¬ 
cused on the performance analysis of wireless systems over 
composite fading channels Q, 0, 0. On a different note, 
diversity techniques are widely used to reduce the effects of 
multipath fading. If the antennas are sufficiently separated, it is 
reasonable to assume that signals received by different antennas 
are independent. However, this assumption is intimately crude for 
systems with closely spaced antennas, such as mobile phones. In 
general, correlation between channels results in a degradation of 
the diversity gain. Therefore, it is important to quantify rigorously 
the effects of correlation in real-life practical scenarios. 

Recently, the Gamma-Gamma (FF) distribution has been pro¬ 
posed as a general and mathematically tractable composite fading 
model 0, 0, 0. The FF distribution is equivalent to the 
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squared generalized-iF distribution and includes the squared K 
distribution and double-Rayleigh distribution as special cases. 
In radio frequency (RF) communications, the performance of 
diversity systems over FF fading channels has been analyzed in 
0, 0. In particular, 0 investigated the outage probability (OP) 
of selection combining (SC) receivers, as well as the average bit 
error rate (BER) of maximal ratio combining (MRC) receivers 
operating over exponentially correlated distributed fading chan¬ 
nels. Furthermore, a simple approximation for the distribution of 
the sum of independent distributed FF random variables (RVs) 
was provided in 0. However, all these works are based on the 
assumption of independent or exponentially correlated FF RVs. 
Furthermore, the FF distribution has been used in the performance 
investigation of free-space optical (FSO) links over atmospheric 
turbulence conditions 0-0. These previous results are not 
theoretically attractive when the underlying FSO links are not 
independently distributed. To the best of the authors’ knowledge, 
a statistical analysis of the multivariate FF distribution with 
arbitrary correlation is still not available in the literature. 

Motivated by the above discussion, this paper makes the 

following specific contributions: 

• Using the Green’s matrix approximation noi, we derived 
new analytical expressions for the probability density func¬ 
tion (PDF), cumulative distribution function (CDF) and mo¬ 
ment generation function (MGF) of the joint FF distribution 
with arbitrary correlation. We point out that the presented re¬ 
sults encompass several previously known results as special 
cases (e.g., those in ||2l, 0, ifTTIl ). 

• Efficient approximations to the PDF and CDF of the sum of 
arbitrarily correlated FF RVs are derived for integer values 
of fading severity m. The resulting expressions can be easily 
and quickly evaluated by using a recursive formula. 

> Furthermore, in order to reveal the importance of the 
proposed statistical formulations, we study the perfor¬ 
mance of SC and MRC receivers over arbitrarily correlated 
generalized-iC fading channels. 

> Finally, in the context of FSO communication systems with 
spatial diversity, we analyze the BER performance when 
strong turbulence channels is assumed. The derived expres¬ 
sion extends the results on independent and exponentially 
correlated FF fading channels in 0. 

Notations: We use upper and lower case boldface to denote 
matrices and vectors, respectively. The expectation is given by 
E(.). The norm of a vector is given by || • ||, while the matrix 
determinant by | • |. The Hermitian operation is defined as 
The set of integer numbers is expressed via Z. Finally, let Xij 
denote the (i, j)-th element of the matrix X. 
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II. Multivariate distribution with arbitrary 

CORRELATION 


A. Multivariate Nakagami-m Distribution with Arbitrary Corre¬ 
lation 


Let r be a Nakagami-m RV, whose PDF is given by nni 


fr (t) 


2^2m-l 

r (m) II™ 



r >0, 


( 1 ) 


where r(-) represents the Gamma function ifTSl Eq. (8.310.1)], 
n = E(r^)/m with E(r^) being the average power, and m > 1/2 
denotes the fading severity. Note that ([T]i is another representation 
of the classical formula for the single Nakagami-m PDF. Without 
loss of generality, let u be a Gamma RV which represents 
shadowing and it is assumed in the following that E(m^) = 1. 
Then, the PDF of u can be expressed as Q 

= u,p>0, (2) 

where /3 is a channel parameter related to the effective number of 
discrete scatterers. Let yi, y 2 , • • ■, y 2 m be A^-dimensional column 
vectors, which are independent of each other and normally 
distributed with zero means and arbitrary correlation matrix S 
which is defined as ^ = 1 for i = j, and ^ = pij for 
i ^ j, where 0 < pij < 1 is the correlation coefficient and 
^ ^ i,j < N. Here we assume that all vectors have the same 
correlation matrix S. 

Now, let xi, X 2 ,..., xat be 2m-dimensional column vectors, 
with x„ composed of the n-th components of the y^ and the 
norm of x„ being r„, i.e., = ||x„||. Then, is a Nakagami- 

m RV. We recall that the joint PDF of N Nakagami-m RVs, 
r = [ri, r 2 ,..., Tat], with arbitrary correlation and identically 
distributed is given by ifTOl Eq. (2)] ifT^ Eq. (9)], as follows 


fr [r) = 


^ PM.Nr% 

H^+™-ir(m) ® 


,2 N-1 


n 


b l 1 —m 

n,n+l| '^n 


X e 




2 \Pn, 


n+l\ 


n 


-r„r„+i 


(3) 


where W denotes the inverse of S with elements pij, and Iy{-) 
is the modified Bessel function of the first kind with order v GSl 
Eq. (8.406)]. 


Remark. Note that, although the parameter 2m seems to be 
restricted as a positive integer, (O can be used for any positive 
values of m not less than 1 /2. This is because it satisfies the 
necessary and sufficient conditions to be a joint distribution 
function, as it was emphasized in m- 

The joint PDE expression Q requires the matrix W to have the 
tridiagonal property. However, in the general case, the inverse of 
S is not a tridiagonal matrix. Therefore, we need to approximate 
E with a Green’s matrix c m, which has the form 


Cit^i 

Cid2 . 

Cl'^N 

Clt^2 

( 2^2 

C2'dN 


C2t^Af ■ 

■ ■ Cn^n 


where Q and iji, are two sequences of real numbers. In order 
to determine C, nonlinear methods (e.g., Levenberg Marquardt, 
quasi Newton, and conjugate gradient) may be used to solve the 


equation C = E, i.e., Q-di — Ei_i = 1 and Qdj = E^j = pij. 
The accuracy of the Green’s matrix approximation technique will 
be verified in Section [V] In any case, its computational com¬ 
plexity, determined by the dimension of the correlation matrix, is 
given by 0{N'^). 


B. Multivariate Gamma Distribution with Arbitrary Correlation 


It is well known that a Nakagami-m RV is the square root 
of a Gamma RV. Furthermore, it has been numerically validated 
in m that the correlation coefficients of Gamma RVs can be 
taken identical to the corresponding coefficients of Nakagami- 
m RVs with satisfactory accuracy. Performing N RV trans¬ 
formations onto ©, the joint PDF of multiple Gamma RVs, 
u = [u;i,aj 2 , ■ • ■, Wat] with uJt = rf = ||xi|p for z = 1,...,iV, 
can be expressed as 


\wre 


- E 


^ (w) = 


—Nm—2 ^ ij 

X Cl 


F (m) 

N-l N-l 

nn 

n=lj=2 




m+zi-1 m+iAT-i —1 




N 


yn-Hj-rHj—i- 




N-l 

n 

n=l 


\Pn, 


n+l\ 


Z„!r (m -I- in) 


■ (5) 


To derive ®, we have used the infinite series representations 
of /«(•) ifTSi Eq. (8.445.1)]. For the special case of exponential 
correlation, where the correlation matrix is defined as E^^ = 1 
for i = j and E^j = pb-il fo]- i f j, ^ reduces to l|5] Eq. (2)]. 


C. Multivariate FF Distribution with Arbitrary Correlation 

The multivariate FF distribution can be derived from the 
product of two multivariate Gamma RVs. Let us assume that 
the Gamma RVs in one multivariate Gamma distribution are 
correlated but they are independent from the Gamma RVs of the 
other multivariate Gamma distribution^ Moreover, let u and lo 
be two independent Gamma RVs with PDFs /«(«) and fui{ijj), 
respectively. The FF RV is denoted by z = uuj. 


Theorem 1. The joint PDF of the multivariate FF distribution 
with arbitrary correlation is given by 


fz {z) = 


X Ff2r,, 


2^\Wf^ 

[r(^)]^r(m) 


E 

n ,...,iAr-i=0 





+ E ^ 




j=l ^3.3 



z„!r(m-|-z„) 


( 6 ) 


where pj = A -^iih foy j = 

aj = iN-ifor j = N, and aj = ij-i-\-ij for j = 2,3,... ,7V—1. 
Also, Kyf) denotes the modified Bessel function of the second 
kind with order v ifTTI Eg. (8.407.1)]. 


Proof: The multivariate FF distribution can be obtained as 

j•ao poo ft / \ 

fz (^)=• y^ n l h dw, n) 

where z = [zi, Z 2 ,..., zn], with z„ = UnU3n- Moreover, (•) is 
the PDF of Un following a Gamma distribution with parameter {3 


'This assumption is reasonable, since in most practical wireless channels the 
effects of fast fading and shadowing are usually independent from each other. 
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given in (2), while u = [ui,U 2 , ■ ■ ■ ,un] is a vector containing N 
independent and identically distributed (i.i.d.) Gamma-distributed 
RVs with the same parameter /3. By substituting (|2]l and (|5]l into 
(|7]i and using ifTSl Eq. (3.471.9)], (|6]l can be derived. ■ 

To the best of the authors’ knowledge, (|6| is new and has not 
been presented in previous literature. Although (|6ll is given in 
terms of infinite series, only a finite number of terms is needed 
to get a satisfactory accuracy (e.g., smaller than 10“®) for all 
considered cases (a detailed discussion is provided in Section FVTi. 


Corollary 1. The joint CDF of the multivariate FT distribution 
z with arbitrary correlation is given by 


F. (z) = 




oo 


[r(/ 3 )]^r(m) 


N 




’ +7 J- 

1 

AT-l 

n 



m + aj, /3,0 

11 

n—1 

i„!F {m+i„) 


where G[-] is a Meijer’s-G function 4751 Eq. (9.301)]. 


( 8 ) 


Proof: We start with the definition of the CDF as 


F. (z) = 


N-1 


xn 


2^\wr 

[r(/3)]'^r(m) 


i„!r (m + in) 


E 

,2Ar_i=0 



2 


N-1 


+ E 

j=i 





f-ij 

Z 

3 

'tj 

Pjj 


K. 


2rij 



dz. 


By transforming the Bessel function Ky (•) into a Meijer’s-G 
function through lIT^ Eq. (14)] and using iflhl Eq. (26)], we can 
conclude the proof with |[^ Eq. (9.31.5)]. ■ 


Note that for (m + Uj — j3) ^ Z and by using ifTTI Eq. 
(07.34.26.0004.01)], ([8]) can be written in terms of the more 
familiar generalized hypergeometric function i7^2 (•) as 


F. (z) = 

N 

n —r 

P 3 ,. 

3=1 


\w\ 


N-1 


- E TT 

( yiP-m-af) ( ppjj 


\Pn, 


ra+l| 


|2i„ 


m+aj 


n 


7„!r (,m+in) 


m+aj 


i^J) ■ ■ 

X 1 F 2 ( m + aj] I — j3+m+aj, l+m+ay, Zj 


r {m+aj— fi) 


1 F 2 I f3', l — m—aj+f3, 1+/3; 


PP3,3 


n 


Corollary 2. The joint MGF of the FF distribution z with 
arbitrary correlation matrix is given by 

\wr 


M, (a) = 


xGli 


N 

[F(/3)]"^F(m) 


— m—a-j 


p, 


’ PP3,3 

1 

iV-1 

n 

' bn.n+ll"*" ' 

Clsj 

m + aj, [3 

11 

n=l 

in!F (m -1- in) 


(9) 


where s = [si, S 2 , • ■ •, sn]- 


Proof: Based on the PDF expression (|6l), the joint MGF of 


the multivariate FF distribution is given by 


nOO nOO j ^ \ 

M^{s)=J ...J expi-'^^SnZn\fz{z)dz. (10) 


Following a similar procedure as in deriving ([8]l, we simply 
express the exponential function in terms of a Meijer’s-G function 
ifThi Eq. (11)] and utilize the relation ifT^ Eq. (21)]. Then, the 
joint MGF expression in (|9|l is derived. ■ 


D. Special Cases 

Now we present some PDF and CDF expressions for special 
cases, e.g., bivariate arbitrarily correlated or exponentially corre¬ 
lated FF RVs. Note that these expressions give the link to previous 
results and prove the correctness of our generalized formulas. For 
bivariate correlated FF RVs, where z = [zi,Z 2 ], the joint PDF 
and CDF reduce to ||2l Eq. (4)] and ^ Eq. (10)], respectively. 
For the special case of the exponentially correlated multivariate 
FF distribution, (|6]l simplifies to ||5] Eq. (5)]. Furthermore, when 
we consider independent FF RVs, the joint PDF, CDF and 
MGF simplify to the product of multivariate independent FF 
distributions ifTTl Eqs. (2), (3) and (4)], respectively. We finally 
point out that for exponentially correlated multivariate FF RVs, 
the joint CDF formula (l8]l reduces to ||5] Eq. (7)] and the joint 
MGF expression ® reduces to ||5] Eq. (10)]. 


III. Sum of FF RVs with Arbitrary Correlation 
A. PDF and CDF of the Sum of FF RVs 
The sum of N FF RVs is defined as 

N N 

S ='^Zn = '^UnUJn, (H) 

n=l n—1 

where are i.i.d. Gamma RVs with parameters {f3, 1//3) whose 
PDF expression was given in (|2|i, and ujn are non-identical and 
correlated Gamma RVs with parameters (to„, n„) and joint PDF 
given in ifTSi Eq. (11)]. According to the approach presented in 
li6l , we can rewrite (dUl as 

^ N N ^ N-1 N 

5'=—^ E (12) 

n—1 n—1 i—1 j—i-\-l 

' ---’ '-V-' 

S e 


where S is the approximation of the sum and e is the approxi¬ 
mation error, respectively. The validity of this approximation has 
been investigated by Kolmogorov-Smirnov (KS) goodness-of-fit 

statistical test in i). Then, S can be expressed as the product of 

N 

two sums of Gamma RVs, i.e. S = S'iS' 2 , where Si = ^ Un 


N 

and 5*2 = X] known that the sum of i.i.d. Gamma 

n—1 

RVs Un remains Gamma distributed with parameters (W/3,1//3) 
lfT9] Eq. (7)]. After variate transformation, we can conjecture that 
Si is Gamma distributed with parameters {Np, 1/N(3). 

Following the similar method presented in ifTSll . let Sn be the 
set of N distinct eigenvalues of the N x N covariance matrix 
given by Ky = E(yyl), where 


F{yt,kyj,i) 


a/2, 

Pi,j /2^ 

0 , 


if i=j and k = l 
if ifj and k = l 
otherwise. 
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We assume that 6n has algebraic multiplicity i/n, where r'„/2 S 
Z. Then, we can derive a closed-form PDF expression for S 2 as 

N rrii 

fS2 (z) = {iJ,{fnn}n=i,{^n}n=l)fu{z-J,ni) , 

i=l j=l 

where /„ (z; j, 17^) is given in (|5]l, m„ = ^'n/2 should be integers, 
= 4:S„/iyn, and 


N 


1 , m / 1 1 ^ 


k ^ ni 

n^i 


X (i,mi - A: +j, 

with /c = 1,..., — 1 and 

N 

n 


(13) 


4'iv (*i trii, , {17^}n=i) — 


n: 




1 1 


With the PDF expressions of S\ and 59, we can derive the PDF 
of 5 all 


N rui 


fs (*’■?’ {mn}n=lA^n}n=l)fl (z; Nj3,j,j^i) 

(14) 

where f^{-) is the PDF of a FF RV and is given by ifTTl Eq. (2)] 
A (z;iV/3,j,ja) = 2 


r(j)r(7V/3) (f7,) 

Finally, the corresponding CDF is given by 


Nl3z\ 

~n~ J 
(15) 


N rui 


“ X/ X/ {'>TT'n}n=lA^n}n=l) p p 


X Gi:i 


'Np 

1 

. 07 ^ 

NP,j,0_ 


(16) 


Note that the approximating PDF and CDF expressions (fl4li and 
(fTST l are both in analytical form, since they involve hnite nested 
summations. The weights coefficients can be easily and quickly 
evaluated by using the recursive formula (fTSl l. 


B. Approximation Accuracy 

In order to evaluate the accuracy of our approximating result, 
the simulated and analytical (O CDF curves of the sum of 
exponentially correlated FF RVs are plotted in Fig. [T] Note 
that the analytical CDF expression can be only obtained by 
calculating the eigenvalues and algebraic multiplicity for each 
specihc covariance matrix. From Fig. [T] it is clear that the 
analytical results approximate the exact ones with good accuracy. 
Moreover, note that the accuracy of the analytical CDF expression 
depends on the values of parameters p, N, /3, m„, and i7„. As a 
general comment, higher values of N and smaller values of i7„ 
yield better accuracy. 

Furthermore, we use a KS goodness-of-ht statistical test to 
validate the accuracy of the approximation results. The KS test 

^ Since the PDF expression involves the distinct eigenvalues of the covariance 
matrix Ky, it is not straightforward to derive similar expressions for the case of 
exponentially correlated FT RVs. Yet, the authors in O have given the PDF and 
CDF expressions for the case of independent FF RVs. 



Fig. 1. Analytical and simulated CDF of the sum of FF RVs (p = 0.1 and 
/3 = 1). 

statistic T is dehned as the maximal difference between the 
simulated and approximate CDFs. The critical value Tmax is given 

’ as Tmax ~ \J §’ where a = 5% is the signihcance level 
and V = 10"* is the number of random samples of RVs ll20ll . Then, 
Tmax = 0.0136. The hypothesis Hq is not always accepted with 
95% signihcance and the accuracy of the proposed approximation 
depends on the combinations of the parameters p, rrin, ^n, 
/3 and N. For example, the hypothesis Hq is accepted with 
T = 0.0118 < Tmax when considering exponential correlation 
with the following typical values, mi = m 2 = 1, f7i = 172 = 1. 
(3 = 1, p = 0.5, and N = 2. However, when we set p = 0.1 and 
keep other parameters the same, the hypothesis Hq is rejected 
with T = 0.0198 > Tmax- 

IV. Applications in Wireless Communications 

A. Diversity Receivers 

In this section, the performance of various classes of multi¬ 
channel diversity receivers, operating over arbitrarily correlated 
generalized-TT fading channels, is analyzed in terms of OP and 
BER. Eor each receive antenna, the transmitted signal is passing 
through the fading channel and perturbed by complex additive 
white Gaussian noise (AWGN) with zero mean and variance Nq. 
The instantaneous signal-to-noise ratio (SNR) of the baseband 
received symbol in the nth diversity branch is given by A„ = 
\hn\^Es/No, where Eg is the energy of the transmitted complex 
symbol and hn is the complex gain of the nth generalized- 
K fading channel. By assuming that each fading channel has 
identical parameters, the average SNR per branch is the same 
and can be expressed as A„ = E^|/i„|^^ Es/Nq = A, where 

E = ml7. Eurthermore, since shadowing occurs in large 

geographical areas, we assume that the shadowing parameters on 
each antenna are same, i.e., = /3, and the AWGN is uncorre¬ 

lated among the diversity channels ED. With A„ = |h„|^A/ml7, 
the joint PDE and CDE of A = [Ai, A 2 ,..., Aat] can be ob¬ 
tained as fxiX) = M^Xi,^X 2 ...,^XN)/{X/mnf 
and Ex (A) = Fz (to 17^, ml7^,..., ml7^), respectively. 

1) Selection Combining: The output SNR of SC diversity 
receivers is the highest instantaneous SNR among the multiple 
branches, i.e., Asc = max{Ai, A 2 ,..., Aat}. An outage event 
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happens when the SNRs of all branches fall below a given 
threshold A*- From (l8]l, the OP is obtained as 


Fout('^th) — 




[r(^)]^r(m) 

Pmpj jXth 




N-1 

n 

' \Pn,r^iy ' 

11 

0 

i„!r {m+in) 


A 


1 

m+aj,l3 ,0 


N 

n —m — OL-j 

j=i 




(17) 


The OP for small values of Ath is typically interesting in 
communication systems. We now focus on the high-SNR regime. 
Based on US] Eq. (9.303)] and lUT] Eq. (07.22.02.0001.01)], the 
high-SNR outage can be derived from (fTTl l and for {m+aj —/3) ^ 
Z as 




IW^I' 


N 


[r(/3)]^r(m; 

N-1 


E n 

j = l 


Tj) 


3,3 




\Pn, 


rJ+l| 


|2i„ 


i„!r (m+i„) 


N 

+0 I I , (18) 


where Tj = min(m + aj,j3) and tj = max(TO + aj,/3). (fTsT l 
indicates that using more antennas at the receiver affects the high- 
SNR slope, since the high-SNR slope grows linearly with N. 

2) Maximal-Ratio Combining: The MRC receiver can achieve 
optimal performance in terms of maximizing the SNR at the 
output of the combiner ||T]. The instantaneous SNR per symbol of 
an 7V-branch MRC receiver is given by Amrc = An- Then, 

the approximating expression for the OP of MRC receivers is 
given by 

N rrii 

-Pout (Ath) {i,j, {^^n}(Ll) 

i=l j=l 

X F.,{X,u;N/3,j,jn,), (19) 

where P’-y(-) is given in (fThl l. Eor (NP — j) ^ Z, the OP in the 
high-SNR regime is given by 

N rrii 

(Ath) = EE^" 

i=l j=l 

^ ^{(^ 3 -Cj) (^n/3^^ ' + o (Ath-^^), (20) 

where Cj = min(7V/3,j) and Cj = max{Nl3,j). It is clear that 
the diversity order depends on the values of the parameters N, (3 
and rrii. Eor the case of MRC, rui in (l20l i should be integer. 

We now consider the average BER of the MRC receiver for 
noncoherent binary frequency-shift keying (NBFSK). Based on 
iU Eq. (9.254)], the average BER expression is given as = 
0.5Mam^c ((?), where g = 0.5 for NBESK. Euthermore, the MGE 
of the output SNR can be obtained from (|9]l as (s) = 

Ma(s, s,...,s). 

Eollowing a similar method of El, the approximate BER of 
MRC receivers for binary phase-shift keying (BPSK) modulation 
is given by 

pMRC ^ (1) + . (21) 

The approximate BER (1211 1 provides a very useful metric of 
communication systems. Eurthermore, (i2Ti) can reduce to El Eq. 
(15)] for exponentially correlated multivariate FF RVs. 


For the BER of BPSK in the high-SNR regime, we invoke the 
seminal parametrization in terms of diversity order Gd and coding 
gain Gc- By keeping only the dominant term of (i2Ti) . the BER of 
MRC receivers over arbitrarily correlated generalized-iT fading 
channels is given by 

pMRC.oo ^ ^ ^ ^ ^ (22) 


where 


N 




(23) 


3 = 1 


Ge = 


iw^r 


N 


N - y TTp7^‘’^'r(f,-T,)r(r,) 


3/3p, 


Ti N-1 


n 


\Pn, 


\2in 


n+l| 


inlF (m -I- in) 


(24) 


Based on the aforementioned expressions, we can notice that 
the diversity order does not depend on the correlation between 
fading channels, but depends on the shape parameter of small- 
scale fading and scale parameter of Gamma shadowing. Note 
that, for large values of /3, the diversity order Gd is limited by 
the shape parameter of small-scale fading. 


B. Free-Space Optical Communications 


We consider a diversity ESO system with N apertures and 
one beam, where the underlying sub-channels between pairs of 
transmit-receive apertures are correlated. Eurthermore, we assume 
that the optimal combining (OC)0 scheme is deployed. Let be 
the fading coefficient of the n-th sub-channel, n = 1,2,..., N, 
that follows the FF distribution with arbitrary correlation. The 
instantaneous electrical SNR at the n-th receive aperture can be 
dehned as A„ = (rjPn)'^/N q and the average electrical SNR 
is Xn = \r]E{Pn)\'^/N q, where -q is the optical-to-electrical 
conversion coefficient. Without loss of generality, we normalize 
the average irradiance as E(P„) = 1 and assume that the average 
electrical SNRs at different receive apertures are equal, i.e., 
Xn — X. The variance of the noise in each receiver is No/2N 
and the average BER is defined as ll^ Eq. (17)] 



vOTVo\ 


N 




fp (p) dp, (25) 


where Q(-) denotes the Gaussian (J-function IT] Eq. (4.1)], p = 
IPi,P2,...,Pm1, and fp (p) is the joint PDE of the vector p 
obtained with the aid of (|6]l. 

In strong atmospheric channels, the fading parameters /3 and 
m can be respectively expressed as IE Eqs. (5)-(6)] 

0.49(7^ 

(l-f 0.18d2-h0.56crf/^y^^ 
0.51cr|(n-0.69crf''®) 

(l -f O.Qd? + 0 . 62 d 2 crf 



/3 = 


exp 


exp 



^In FSO communications, the term OC is used instead of MRC. 
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Fig. 2. Analytical, simulated and high-SNR approximation OP of SC diversity 
receivers with three branches over exponentially correlated generalized-iC fading 
channels against the inverse normalized outage threshold A/Ath (p = 0.25 and 
m = 2). 


where cr| = and d = y/kD'^/4L with L being 

the distance between transmitter and receiver, k is the optical 
wavenumber, D is the aperture diameter and = 1.7 x 10“^"* 
is the refractive-index structure parameter. 

The integral in (l25T l is very difficult to be evaluated in closed- 
form, so we use a simple and accurate exponential approximation 
for Q{-) (si Eq. (14)] to derive an approximate result. Then, (l25T l 
can be written as 


P 


oc 


12 


9 N 

V- o2 


ANNo 


n—l 


■ exp - 


ri 


N 


3NNo 


E pp ^p- 


(28) 


Substituting (|6ll into (l28T l and using ifThl Eqs. (11), (14)] and ifThl 
Eq. (21)], we derive the average BER as 




2^ 

4/V 


where 


A (m, /3, x) = 


N 




2Af(m+/3-2)|p^|»^ 

4xfl \ ^ 


m,/3, ^ 

/ A \ 

'mj’ 

N-1 ■ 

TT 

|2/?n,n-|-l 


(29) 


i„!r (m+in) 




1—(rn+ctj) 2-{m+aj) l — j3 2—j9 
2 ’ 2 > 2 > 2 
0 


Note that for the special case of exponentially correlated 
atmospheric turbulence, ( |29] | coincides with a previously known 
expression 0 Eq. (19)]. 


V. Numerical Results 
A. RF Communication Systems 

In our simulation, without loss of generality, it is assumed 
here that Ii = I 2 = ... = In-i and every branch has the 
same average SNR A. Eor the OP of three branches with SC 
receivers, the Monte-Carlo simulation results, analytical results 
dnii and high-SNR approximation (fTsT i are compared in Eig. 
121 We again assume exponential correlation. The shadowing 
parameter f3 corresponds to 1.0931, 7.9115 and 75.1155 for 



Fig. 3. Average BER of triple SC diversity receivers over linearly correlated 
generalized-iC fading channels against the average SNR per branch A. 



Fig. 4. OP of MRC diversity receivers over exponentially correlated generalized- 
K fading channels against the inverse normalized outage threshold A/Ath W — 
ran = m = 1). 

frequent heavy shadowing, Karasawa shadowing and infrequent 
light shadowing, respectively 0. It is clear that for exponentially 
correlated generalized-AT fading channels, the analytical results 
agree closely with the exact OP results obtained via simulations. 
Eurthermore, the high-SNR approximations are sufficiently tight 
and become exact, even at moderate SNR values. As expected, the 
OP that corresponds to light shadowing (/3 = 75.1155) is larger 
than the OP that corresponds to heavy shadowing (/3 = 1.0931). 
Moreover, this trend is more pronounced for smaller values of /3. 

Eigure [3] presents the simulated, analytical (12111 and high- 
SNR asymptotic results (l22ll for the BER of triple SC diversity 
receivers over linearly correlated generalized-Tf fading channels. 
The linearly correlated model is a more general case than the ex¬ 
ponential one with a correlation Toeplitz structure matrix defined 
as = 1 for i = j and = pij for i ^ j. We note 

that as the scale and shadowing coefficients increase, the average 
BER decreases. Eor all scenarios, the analytical results agree with 
the Monte-Carlo ones. Likewise, the diversity order and coding 
gains are accurately predicted. It is clear that the diversity order 
depends on the shadowing parameters. 

The simulated and approximate OP of MRC receivers over 
exponentially correlated generalized-AT fading channels are inves- 
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Fig. 5. Bit eiTor rate of quadruple branch MRC diversity receivers over 
exponentially con'elated generalized-iC' fading channels against the average SNR 
per branch A (m = 2, p = 0.3). 



Fig. 6. Average BER of space diversity FSO system with three apertures and 
one beam over strong turbulence fading channels against the average electrical 
SNR A. 

tigated in Fig. |4] It is obvious that the proposed approximation is 
a lower bound of the simulated results. However, the difference is 
not greater than 2 dB in all cases. Furthermore, the approximation 
results are more accurate with the increased correlation coefficient 
p and the decreased value of N. 

We present the simulated and approximate BER perfor¬ 
mance of quadruple MRC receivers over exponentially correlated 
generalized-^ fading channels in Fig. |5] It is clear that the BER 
decreases as /3 increases. Approximate BER curves are close 
to the simulated results and can be used as an upper bound. 
Moreover, the gap between two curves becomes smaller as the 
SNR increasing. 

B. FSO Systems 

In Eig.|6] the simulated and approximate (|29] | error performance 
of space diversity ESO links with N = 3 receive apertures 
employing OC over exponentially correlated atmospheric turbu¬ 
lence channels, is depicted. We consider the case of the optical 
wavenumber k = 0.405 x 10^, link distances L = 3000, 5000m 


and correlated coefficients p = 0.2,0.5,0.7, respectively. Upon 
using (|2^ and dZTl l. the respective values for /3 and m can be 
determined. Erom Eig. |6l it is clear that the approximate results 
of cases under consideration closely agree with the ones obtained 
via Monte-Carlo simulations. Note that the average BER increases 
by increasing p and L or decreasing A. 
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